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We review the recent glueball mass calculations using an efficient method for solving the Schrodinger equation 
order by order with a scheme preserving the continuum limit. The reliability of the method is further supported 
by new accurate results for (l+l)-dimensional a models and (2+l)-dimensional non-abelian models. We present 
first and encouraging data for the glueball masses in 3+1 dimensional QCD. 



1. INTRODUCTION 

Although lattice field theory is expected to give 
the most reliable estimates for QCD spectroscopy, 
the values for glueball masses have still been an is- 
sue under debate j|] for the last two decades. Re- 
cently, more accurate numerical calculations by 
the IBM group |^ on much larger lattices and 
higher statistics lead to M(0++) « 1.740 Gev in 
the infinite volume and continuum limit. 

Alternatively, the lattice Hamiltonian meth- 
ods jH"^ receive increasing attention, since more 
physical insights can be obtained and significant 
progress has been achieved in recent years. The 
gluon dynamics is described by the Hamiltonian 
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The most direct way to extract both the glue- 
ball masses and their wavefunctions is to solve the 
Schrodinger eigenvalue equation H\F) = epW). 
The glueball wavefunction \F) is created by the 
gluonic operators with the given J^'-' acting on 
the vacuum jfi). The vacuum wavefunction jfJ) 
satisfies H\Q) = en|ri). An estimate of the glue- 



ball mass is then M jrc 
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Physically, the low energy spectrum originates 
mainly from the long wavelength excitations. In 
a series of papers we developed further a 

method for solving the eigenvalue problem in 
a new scheme preserving the correct long wave- 
length behavior at any order of approximation. 
Having such a correct truncation is essential for 



getting a correct continuum behavior of the phys- 
ical quantities. This idea and the scheme are con- 
firmed by the results of three-dimensional abelian 
in 1^ and non-abelian in U-|^Jl^ gauge theories: 
they converge rapidly, and even at very low trun- 
cation orders clear scaling windows for the vac- 
uum wavefunction and mass gaps have been es- 
tablished. For 2+1 D U(l) and 2+1 D SU(2), 
they agree perfectly with recent accurate Monte 
Carlo data. In our pioneering study of 2+1 D 
SU(3), we obtained the first estimates for the vac- 
uum wavefunction |^ and the glueball masses . 
Most of these results and detailed techniques have 
been summarized last year JTT| , p^ . 

There have been new developments this year. 
Increasing evidences, higher order and more ac- 
curate data, including those [|l3| for the two- 
dimensional 0{N) a models, have been obtained 
to support the reliability of the approach. Here 
we would like to highlight our recent results 
||l|,|l|l for QCD in 2+1 and 3+1 dimensions. 



2. QCD3: BEYOND A TOY MODEL 

The reader will amaze later at the fact that 
QCD3 is not just a toy model of QCD4, but it can 
even mimic the realistic theory. Using the tech- 
niques in 16 1 and analysing carefully the data 
in the observed scaling region /9 G [5, 12), a more 
accurate value jlj] for AI{0'^'^)/e'^ is obtained: 



M(0++) 



2.15 ±0.06, 
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with the error being the systematic uncertainties 
due to finite order truncation (in it was 2.1 
for narrower /? range at third order). (H) can be 
compared with Samuel's result 1.84 ± 0.46 in the 
continuum Hamiltonian formulation or the 
preliminary MC result 2.4 ± 0.2 on a finte lattice 
obtained much later 10. 

A relation between the glueball mass and the 
confinement scale from the vacuum wavefunc- 
tion may be induced. The vacuum functional, 
which interpolates the strong and weak coupling 
regimes, are [ p^|Jl7| 

\n) = exp{^ J d^-ix tr[T{V'+rT'^^^]}. 

The correlation length ^, with dimension of in- 
verse mass, is proportional to e~^, i.e., the con- 
finement scale in the vacuum. might also be 
relevant for the constituent gluon mass and the 
lightest glueball mass ||l^. In the strong cou- 
pling or large Nc limit, I f2) reduces to the strong 
coupling wavefunction ||3|,|20|| . In the intermediate 
and weak coupling, it becomes [p^PJEfl 

\n) = exp{ J d^-^x [-^lotrJ'^ - H2 triVJ'f]}, 

identical to our long wavelength vacuum wave- 
function I^J^. The correlation length is then re- 
lated to fio and fj,2 hy ^ — (—2^2/^0)^^^- For 
2+1 D SU(2), ^ = 0.65/e2 (see |l|,|), while 
for 2+1 D SU(3), our result g is ^ = O.SS/e^. 
If the glueball mass is proportional to the con- 
stituent gluon mass, from the difference of the 
scales between SU(2) and SU(3), one may also 
guess M(0++)/e^ « 2, consistent with (||) from 
our practical calculation. 

Combining the recent MC data |2T] for the 
string tension a in QCD3, we obtain in the con- 
tinuum limit 
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3.88 ±0.11. 
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To extend QCD3 to QCD4, we follow the di- 
mensional reduction argument ||Jl^, which says 
that a confining theory in D dimensions {2 < D < 
4) becomes a localized field theory in = I? — 1 
dimensions. This can be exactly proven in the 
strong coupling or large Nc limit. In this limit. 



the fixed time vacuum expectation value of an 
operator 0(U) in D dimensions corresponds to 
the path integral expression for < 0{U) > in 
D — 1 dimensional lattice field theory. In the 
intermediate coupling region the 3+1 D theory 
can still be approximated by its 2+1 D theory for 
long wavelength configurations in comparison to 
the confinement scale. Accordingly, M{J^'~^)/ y/a 
for the lightest glueball should be approximately 
the same for 2+1 and 3+1 dimensions. Since for 
SU(3), Nc is larger and the measured length ^ in 
the vacuum functional (^) is smaller than those 
for SU(2), our speculation is that dimensional re- 
duction work better. In fact, (^ is well consis- 
tent with the IBM data M(0++)/Vct = 3.95 from 
MC simulation of QCD4, providing the contin- 
uum = 0.44 Gev is used. From this world 
averaged ^/a and (|^), we expect 



A-f (0++) = 1.71 ± 0.05 GeV, 



(4) 



in nice agreement with the IBM data M(0++) = 
1.740 ± 0.071 H. This favors 6'//j(1710) as a 
candidate of the 0++ glueball. 

3. QCD4: THE REALISTIC THEORY 

It is very desirable to do concrete computa- 
tions in the realistic theory: QCD in 3+1 dimen- 
sions. In |]l^, we made a first step towards this 
direction. As our first attempt, we computed the 
masses of glueballs 0"''"'", 0~~, and 1"' , which 
gluonic operators are easily constructed. 

At relatively strong coupling, the absolute 
value of uMe, calculated in the Hamiltonian for- 
mulation and converted to the Euclidean one, dif- 
fers from the result from Euclidean strong cou- 
pling expansion. This is not surprising because 
they are different schemes at finite lattice spac- 
ing, and the weak coupling relation between them 
doesn't not hold for strong coupling. For (3 > 6.0, 
we do observe much smaller difference in aMs 
and a tendency approaching the MC data . 

Similar to the most recent MC data [Q in the 
available coupling region, clear asymptotic scal- 
ing window for the individual mass aM{0'^~^), 
aM(0 ), or aM(l+") could not be found. 
There might be two possible reasons for this scal- 
ing violation: the available coupling region be- 
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ing not weak enough for the asymptotic scal- 
ing law to be valid or the results being not ac- 
curate enough. The first one may be reduced 
by the Symanzik's improvement or the Lepage- 
Mackenzie scheme. The second one may be im- 
proved by higher order calculations. For the mass 
ratio M^jp /M''j^\ however, the large part of er- 



rors m m]^^^ and M)^> are cancelled. Indeed, we 
observe approximate plateau in the mass ratio 
both for Af(0~)/M(0++) and M{1+') / M{0++) 
in the coupling region e [6.0, 6.4]. 

From the plateaus in /3 e [6.0, 6.4], we estimate 
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M(l+-) 
M(0++) 



= 2.44 ±0.05 ±0.20, 



= 1.91 ±0.05 ±0.12, 



(5) 



where the mean value is the averaged one over 
the data in this region, the first error is the error 
of the data in the plateau, and the second error 
is the upper limit of the systematic errors due to 
the finite order truncation. 

For comparison, we list the corresponding 
results from other lattice calculations. For 
M(0~)/M(0++), it is approximately 3.43 ±1.50 
from Monte Carlo simulations with a huge error 
bar (see the references in j^]), and 2.2 ± 0.2 from 
t-expansion plus Fade approximants ] p^ . For 
M(l+-)/M(0++), it is 1.88 ± 0.02 from Monte 
Carlo, and 1.60±0.60 from t-expansion plus Fade. 
For M(0 )/M(0++), we obtain a value between 
the Monte Carlo and t-expansion data, with error 
under much better control than the former one. 
For the M(l+-)/Af(0++), where the systematic 
error in the MC data is very small, we get a value 
consistent well with them. 

To summarize, we have tested with higher pre- 
cision our new scheme [Q-^ in 1+1 D a models 
and 2+1 D non-abelian models. For the light- 
est glueball mass in QCD3, accurate results are 
obtained with systematic uncertainty under well 
control. The idea of dimensional reduction is suc- 
cessively applied to extrapolate QCD3 to QCD4. 
We have also done concrete calculations of QCD 
in 3+1 D and presented first and encouraging 
data for some glueball masses. The inclusion of 



other glueballs, such as 2++, 0""*" and O"""", and 
reduction of the systematic errors are in progress. 
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